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ABSTRACT 


An  abstract  approximation  framework  for  the  identification  of 
nonlinear  distributed  parameter  systems  is  developed.  Inverse  problems  for 
nonlinear  systems  governed  by  strongly  maximal  monotone  operators 
(satisfying  a  mild  continuous  dependence  condition  with  respect  to  the 
unknown  parameters  to  be  identified)  are  treated.  Convergence  of 
Galerkin  approximations  and  the  corresponding  solutions  of  finite 
dimensional  approximating  identification  problems  to  a  solution  of  the 
original  infinite  dimensional  identification  problem  is  demonstrated  using 
the  theory  of  nonlinear  evolution  systems  and  a  nonlinear  analog  of  the 
Trotter-Kato  approximation  result  for  semigroups  of  bounded  linear 
operators.  The  nonlinear  theory  developed  here  is  shown  to  subsume  an 
existing  linear  theory  as  a  special  case.  It  is  also  shown  to  be  applicable 
to  a  broad  class  of  nonlinear  elliptic  operators  and  the  corresponding 
nonlinear  parabolic  partial  differential  equations  to  which  they  lead.  An 
application  of  the  theory  to  a  quasilinear  model  for  heat  conduction  or 
mass  transfer  is  discussed. 
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1.  Introduction 

In  this  paper  we  develop  a  general  abstract  approximation  framework 
for  the  identification  of  nonlinear  distributed  parameter  evolution  systems. 
Our  intent  is  to  define  relatively  straightforward  and  easily  verified 
criteria  that  are  applicable  to  broad  classes  of  nonlinear  systems;  these 
criteria  will  guarantee  the  convergence  of  solutions  to  a  sequence  of  finite 
dimensional  Galerkin  approximation  based  parameter  estimation  problems  to 
a  solution  of  the  original,  underlying,  infinite  dimensional  identification 

f 

problem.  The  results  that  we  present  below  generalize  and  extend  the 
theory  recently  developed  by  Banks  and  Ito  in  [2]  and  [3]  for  regularly 
dissipative  or  abstract  parabolic,  linear  systems.  It  is,  to  the  best  of  our 
knowledge,  the  first  such  general  approximation  theory  for  inverse 
problems  involving  nonlinear  distributed  systems. 

The  sufficient  conditions  set  down  in  our  framework  include  a 
relatively  mild  continuity  assumption  with  respect  to  the  unknown 
parameters  to  be  identified,  an  equi-boundedness  and  an  equi-strong 
monotonicity  assumption  on  the  nonlinear  operator  describing  the  system 
dynamics.  In  addition  our  theory  requires  a  standard  approximation 
assumption  on  the  Galerkin  subspaces  used  to  effect  the  finite  dimensional, 
or  finite  element,  approximations.  We  demonstrate  that  solutions  to  the 
finite  dimensional  identification  problems  approximate  a  solution  to  the 
infinite  dimensional  identification  problem  via  a  convergence  result  for 
solutions  to  the  forward  problems.  This  result  is  obtained  using  the  theory 
of  nonlinear  evolution  systems  and  a  nonlinear  analog  of  the  well-known 
Trotter-  Kato  approximation  result  for  linear  semigroups. 


In  the  present  paper,  we  are  concerned  only  with  theory; 
implementation  questions  and  conclusions  drawn  from  our  numerical  or 
computational  studies  will  be  reported  on  elsewhere.  Also,  while  we  have 
tried  to  make  our  framework  as  versitile  as  possible,  the  treatment  below 
does  have  limitations.  For  example,  our  theory  can  handle  quasi-autonomous 
systems  but  it  is  not  applicable  in  the  fully  nonautonomous  case.  The 
development  of  a  general  theory  which  can  handle  nonlinear  systems 

involving  time  dependent  operators  requires  additional  effort  and  is 
currently  the  focus  of  our  ongoing  investigations.  The  particular 
difficulties  inherent  in  the  time  dependent  case  will  be  described  in 
greater  detail  in  our  discussions  below. 

We  provide  a  brief  outline  of  the  remainder  of  the  paper.  In  Section 
2  we  state  a  fundamental  existence  and  uniqueness  result  for  infinite 

dimensional  nonlinear  systems  and  prove  a  general  approximation  result 
which  is  especially  well  suited  for  application  in  the  context  of  the  inverse 
problems  which  are  the  central  focus  of  our  study.  In  Section  3  we 

define  a  class  of  nonlinear  distributed  systems  and  the  associated 
parameter  identification  problems.  We  define  the  Galerkin  approximations 
and  prove  the  general  convergence  result.  Section  4  contains  some  examples. 
We  show  that  our  nonlinear  theory  subsumes  the  linear  theory  presented  in 
[2J  and  [3]  as  a  special  case;  we  also  consider  the  application  of  our 

framework  to  a  class  of  nonlinear  elliptic  operators  and  the  corresponding 
nonlinear  parabolic  partial  differential  equations  to  which  they  lead.  In 
particular,  we  look  at  the  application  of  our  results  to  a  well  known 
quasilinear  model  for  heat  conduction  or  mass  transfer.  In  Section  S  we 
summarize  our  findings  and  provide  some  concluding  remarks. 
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2.  An  Approximation  Result  for  Nonlinear  Evolution  Systems 

Let  X0  be  a  Banach  space  with  norm  |-|0.  We  consider  the  nonlinear, 
quasiautonomous  initial  value  problem  in  X0  given  by 

(2.1)  x0(t)  +  A0x0(t)  *  f0(t),  0  <  t  <  T, 

(2.2)  x0(0)  -  x° 

x 

where  x®  €  X0,  f0:  [0,T]  -»  X0  and  the  nonlinear  operator  A<y  X0  ■*  2  0  is  in 
general  multivalued,  not  everywhere  defined,  and  not  continuous.  The 
existence  of  solutions  to  the  initial  value  problem  (2.1),  (2.2)  and  the 
subsequent  approximation  result  to  follow,  are  both  consequences  of  Theorem 
2.1  to  be  given  below. 

We  shall  require  the  following  definitions.  Let  X  be  a  Banach  space  with 
norm  |-|x.  For  A:  X  -*  2X,  a  nonlinear,  multivalued  operator,  the  domain  and 
range  of  A  are  defined  by  Dom(A)  -  {x  €  X:  Ax  *  0}  and  R(A)  «  u  Ax 

X€dom(A) 

respectively.  We  say  that  the  operator  A  is  accretive  if  for  every  X  >  0, 
Xj.Xj  €  Dom(A)  and  yj  €  Axr  ya  €  Axa  we  have 


I 


IXj  -  X,|x  <  |Xj  -  Xa  +  V(ya  -  ya)|x. 

We  say  that  A  is  m-accretive  if  A  is  accretive  and  X(I  +  XA)  «  X  for  some 
X  >  0.  We  note  that  if  A  is  m-accretive  then  R(I  +  XA)  >*  X  for  every  X  >  0 
and  for  each  X  >  0  the  resolvent  of  A  at  X,  J(X;A):  X  -»  X,  a  single  valued, 
everywhere  defined,  nonlinear  operator  on  X  can  be  defined  as  J(X;A)  * 
(I  +  XA)'1. 


A  two  parameter  family  of  nonlinear  operators  (U(t,s):  0  <  s  <  t  <  T) 
defined  on  a  subset  Q  C  X  is  called  a  nonlinear  evolution  system  on  ft  if  for 
each  x  €  ft  we  have  U(t,s)x  €  ft,  U(s,s)x  *  x  and  U(t,r)U(r,s)x  «  U(t,s)x  for 
0<s<r<t<T  and  U(t,s)x  is  continuous  from  the  triangle  A  =  {[s,t]: 
0  <  s  <  t  <  T}  into  X. 

A  strongly  continuous  function  x:  [0,t]  -  X  is  called  a  strong  solution  to 
the  quasiautonomous  initial  value  problem 

(2.3)  x(t)  +  Ax(t)  3  f(t),  0  <  t  <  T 

(2.4)  x(0)  -  x° 

where  f:  [0,T]  -•  X  and  x°  €  X  if  x  is  absolutely  continuous  on  compact 
subintervals  of  (0,T),  differentiable  almost  everywhere  and  satisfies  f(t)  - 
x(t)  €  Ax(t)  for  almost  every  t  €  [0,T]  and  x(0)  *  x°. 

Theorem  2.1.  Let  X  be  a  Banach  space  with  norm  |-|x  and  suppose  that  A:  X  -* 
2X  and  f:  [0,T]  -•  X  appearing  in  (2.3)  satisfy 

(1)  there  exists  anu€  R  for  which  the  operator  A  +  u>I  is  m-accretive, 

(2)  feL^O.TiX). 

Then  a  unique,  nonlinear  evolution  system  (U(t,s):  0  <  s  <  t  <  T }  on  Dom(A)  can 
be  constructed  which  satisfies 

(i)  |U(t,s)$  -  U(t,s)<flx  <  eu(t*,)|0  -  #x,  for  $,il>  €  Dom(A)  and 
0  <  s  <  t  <  T, 

(ii)  |U(s+t,s)$  -  U(r+t,r)<«x  <  2  [  ew<*-T>|f(T+s)  -  ^T+rJI^T, 
for  all  $  €  Dom(A)  and  all  t  >  0  such  that  s+t,  r+t<T. 
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(iii)  if  x°  €  Dom(  A)  and  the  initial  value  problem  (2.3),  (2.4)  has  a  strong 
solution  x,  then 

x(t)  =  U(t,s)x(s),  for  0  <  s  <  t  <  T. 

When  x°  €  Dom(A),  the  strongly  continuous  function  x:  [0,T]  -*  X  given  by  x(t) 
U(t,0)x°  is  referred  to  as  a  mild  or  generalized  solution  to  (2.3),  (2.4). 

Theorem  2.1  is  a  direct  consequence  of  results  given  by  Crandall  and 
Evans  and  Evans  in  [7]  and  [9].  Henceforth,  we  shall  assume  that  A0:  XQ  -»  2xo 
and  f0:  [0,T]  -•  X,  satisfy  (1)  and  (2)  in  the  statement  of  Theorem  2.1  and 

that  x0  €  Dom(A0).  We  then  let  (U0(t,s):  0  <  s  <  t  <  T)  denote  the 
corresponding  nonlinear  evolution  system  on  Dom(AQ)  and  consider  the 
approximation  of  mild  solutions  to  the  initial  value  problem  (2.1),  (2.2). 

Our  approximation  result  is  in  the  spirit  of  those  given  for  nonlinear 
semigroups  and  evolution  systems  by  Crandall  and  Pazy  in  [8]  and  Goldstein 
in  {10],  However,  our  theorem  differs  from  these  earlier  treatments  in  two 
ways.  First,  we  require  that  the  time  dependent  perturbation  f0  be  only  Lt  as 
opposed  to  it  being  continuous  as  in  [8]  and  it  satisfying  a  Lipschitz-like 
condition  in  [10].  This  distinction  is  especially  relevant  in  the  case  of  control 
systems  where  discontinuous  input  is  common.  The  second  difference  is  that 
we  give  our  result  in  a  form  that  is  most  appropriate  for  application  to  the 
development  of  a  general  approximation  theory  or  framework  and 
computational  schemes  for  the  parameter  identification  problems  to  be 
discussed  in  the  next  section. 

We  shall  require  some  set  theoretic  notation.  For  sets  Hn,  n»0,l,2,...,  by 
lim  HnDH0we  shall  mean:  Given  x0€H0,  there  exist  xn€Hnsuch  that  xn-x0as 
n  ->  •. 


Theorem  2.2.  For  each  n  €  Z+  *  {1,2,3, ...)  let  Xn  be  a  closed  linear  subspace  of  XQ. 

x 

For  n-  0,1,...,  let  An:  Xn  -*  2  n  be  a  possibly  multivalued  nonlinear  operator  on  Xn,and 
let  fn:[0,T  ]-*  Xn  bean  Xn-valued  measurable  function  defined  on  [0,TJ.  Suppose 
that  there  exists  an  u0  €  R,  independent  of  n,  for  which  the  operators  An  +  uQI 
are  m-accretive ,  that  there  existsa  function  g€Lt  (0,T;5^)  for  which  |f n (t)|  <  g(t), 
a.e.t€  [0  ,T],  and  that  lim  Dn  D  DQ  where  Dn  =  Dom(An)  and  D0  =  Dom(A0). 
Suppose  further  that  for  some  XQ  >  0  we  have 

(2.5)  lira  J(XQ;  A„  +  wQl)*D  =  J(XQ;  AQ  +  u>Ql)4>0 

whenever  4>  €  X  with  lim  6  ■  €  Xft  , 

n  n  n_M»  n  0  O’ 

and  that 

lim  fn(t)  =  f0(t)  for  a.e.t€  [0  ,TJ. 

Then  for  each  n  €  Z+  there  exists  a  unique  nonlinear  evolution  system  {Un(t,s): 

0  <  s  <  t  <  T}  on  Dn  corresponding  (in  the  sense  of  Theorem  2.1)  to  An  and  fn 

end  for  <bn  €  Dn  with  lim  4>n  m  $0  G  DQ  we  have 
n"4® 

(2.6)  Hm  U„(t,s)^n  =  Uo(t,s)0o,  0  <  s  <  t  <  T, 
with  the  limit  being  uniform  in  t  for  t  €  [s,T], 

Proof.  We  follow  Goldstein  (see  [10],  [11])  and  use  an  approach  first  suggested 
by  Kisynski  [13]  for  demonstrating  the  convergence  of  approximations  to 
linear  semigroups,  to  prove  the  theorem  via  an  application  of  our  existence 
result.  Theorem  2.1. 
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Let  X  =  {x  *  {x„C=0:  Xn  6  Xn*  n  “  0*1>2 .  and  Xn  =  X0>  and  for 

x  €  x  set  Ixl  =  sup|xn|0.  Then  l-l  defines  a  norm  on  the  linear  vector  space 

n 

X,  and  the  space  x  together  with  the  norm  1-1  is  a  Banach  space.  Define  the 
operator  A:  X  •*  2X  by 

dom(^)  =  {x  =  {xn}“_0  €  X:  xn  €  Dom(An),  and  for  each  n  «=  1,2,  ... 
there  exists  a  yn  €  Anxn  such  that  lim  yn  «  y0  €  A0x0), 

n"*® 

for  x  €  Dom(/4),  y  =  {yn}“=o  €  ^x  if  and  only  if  y„  €  A„xn’ 
n  =  0,1,2,  ...  and  lim  yn  =  yQ. 

n*4® 


Define  an  essentially  x-valued  function  f  on  the  interval  [0,T]  by  f(t)  = 
{fn(t)}*_o-  The  assumptions  on  the  fn  are  such  that  fn(t)  -*  f0(t)  for  almost 
every  t  €  [0,T],  However,  by  appropriately  redefining  on  a  set  of  measure 
zero,  we  may  infer  from  the  assumptions  on  the  functions  fn  that  f:[0,T]  -*  x 
with  f  €  LjfO.TiX). 

It  is  readily  seen  that  the  operator  A  +  u>0I  is  m-accretive.  Let  x1  = 
{xi>"=0,  x2  =  {**}"=0  €  Dom(/l)  and  let  y1  -  {yx}~=0  e  Ax1  and  y2  -  {y2}~=0  € 

A  « 

Ax  .  Since  for  each  n  =  0,1,2,  ...,  An  +  u0I  is  assumed  to  be  m-accretive,  for 
X  >  0  we  have 


lx1  -  x2l  =  sup  |x*  -  x2|0  <  sup|x*  -  x2  +  X(yi 


n  +  U0X„ 


1  _ 


Cy  l 


U0Xn 


*))lo 


lx 


1  a2 
1  -  X*  + 


My1  +  WqX1  -  (y* 


uox 


*))■ 


and  therefore  that  A  +  wQI  is  accretive.  Now  let  y  «  (yn}"_0  €  x  and  sct  x  * 

<XX=0with  Xn  *  J<*0;  An  +  ^o^n’  n  “  ~  whcre  X0  is  chosen  as  in  (2-5).  It 

is  immediately  clear  that  for  each  n  ■=  0,1,2,  ..., 


xn  €  Dom(An)  C  Xn.  Since 


y  €  x  we  have  lim  y  »  y0  and  therefore,  by  assumption  (2.5),  that  iim  x  =  x. 
or  x  €  x  Setting  zn  -  (yn  -  (1  +  n  -  0,1,2,  ....  it  follows  that 

z  €  A  x  and  lim  z_  -  zA  €  A„x„.  We  conclude  that  x  €  Dom(/4),  (I  + 

n  n  n  °  °  u  u 

x0(^+u0I))x  a  y,  and  that  X(I  +  x0(^  +  w0I))  -  x 

We  have  shown  that  the  operator  A  and  the  function  f  satisfy  conditions 
(1)  and  (2)  given  in  the  statement  of  Theorem  2.1.  Therefore,  a  unique 
nonlinear  evolution  system  (U(t,s):  0  <  s  <  t  <  T)  on  Dom(^)  corresponding  to 
A  and  f  can  be  constructed  with  U(t, s)  -  (t/n(t,s)}"_o-  Using  assumption  (2.5) 
it  can  be  shown  that  Dom(A)  -  (x  -  (xn}"=o  €  )C  xn  €  Dn,  n  «=  0,1,2,  ...  and 
lim  x  -  xj.  Since  X(U( t,s))  C  x  it  follows  that 

(2.7)  lim  f/n(t,s)#n  -  CUMXL  0  <  s  <  t  <  T 
n'4m 

whenever  $n  €  Dn  and  lim  ■  $0  €  D^.  Since  each  of  the  operators  An  and 

n 

the  functions  fB  satisfy  conditions  (I)  and  (2)  of  Theorem  2.1,  unique 

nonlinear  evolution  systems  (Un(t,s):  0  <  s  <  t  <  T)  on  Dn  corresponding  to  An 

«  _ 

and  f  can  be  constructed.  Recalling  that  Dom(>4)  C  x  D  ,  we  may 

“  n=0  “ 

define  the  family  of  operators  (V(t,s):  0  <  s  <  t  <  T)  on  Dom(  A)  by 

(2.8)  V(t,s)J  -  <Vn(t,s)xn>-=0  s  (UJt.sJxJ-^ 

for  x  -  {Xn)“=0  €  Dom(/4).  Uniqueness  (see  [9])  dictates  that  for  each  n  « 

0,1,2 .  C/„(t,s)xn  -  Vn(t,s)xn  whenever  (xn}"=0  €  Dom(^).  This  together 

with  (2.7)  and  (2.8)  establish  (2.6).  The  fact  that  the  convergence  in  (2.6) 
is  uniform  in  t  for  t  €  (s,T)  is  argued  exactly  as  it  was  for  the 
convergence  of  approximations  to  nonlinear  semigroups  in  the  proof  of 
Theorem  3.2  in  [10]. 
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We  note  that  (2.S)  is  also  a  necessary  condition  for  the  conclusion  to 
hold  (see,  for  example.  Theorem  1  in  [14]). 


3.  An  Approximation  Theory  for  Identification  Problems 

Let  H  be  a  real  Hilbert  space  with  inner  product  <•,•>  and 

corresponding  norm  |  |.  Let  V  be  a  reflexive  real  Banach  space  with 

norm  I  •  I  and  let  V*  be  its  dual.  (All  of  our  theory  can  be  developed  in 

complex  spaces  if  necessary;  see  [6].)  We  denote  the  usual  dual  norm  on  V* 

by  I -I*  and  assume  that  V  is  densely  and  continuously  embedded  in  H 
with  |v|  <  4I  v  I,  v  €  V,  for  some  positive  constant  n.  Identifying  H  with 
its  dual,  we  obtain  V  c  H  -  H*  c  V*.  For  4>  €  V*  and  v  €  V  the  duality 
pairing  between  $  and  v  is  denoted  by  <$,v>.  When  $  €  H,  its  pairing 
with  v  €  V  agrees  with  the  inner  product  of  $  with  v.  It  follows  for  u  e 
H  and  v  €  V  that  lul*  <  *t|u|  and  Ivl,  <  I v  1 .  Let  Q  and  Z  be  metric 
spaces  and  let  Q  be  a  nonempty,  sequentially  compact,  subset  of  Q.  The 
spaces  Q  and  Z,  and  the  set  Q  are  referred  to  as  the  parameter  space,  the 
observation  space,  and  the  admissible  parameter  set  respectively. 

We  recall  that  a  single  valued  operator  A:V-.V*  is  hemicontinuous  if 

lim  A(u+tv)  ■  Au  for  all  u,v€V  where  the  limit  is  taken  in  the  weak  sense. 
t-«o 

For  each  q  €  Q  let  A(q):  V  -»  V*  be  a  single  valued,  hemicontinuous,  (in 
general,  nonlinear)  operator  satisfying: 

(A)  (Continuity):  For  each  v  €  V,  the  map  q  -*  A(q)v  is  continuous  from 
Q  C  Q  into  V*. 

(B)  (Equi  V-  monotonicity):  There  exist  an  w  €  R  and  an  «  >  0,  both 
independent  of  q  €  Q,  such  that 

<A(q)u  -  A(q)v,  u-v>  +  utu-v|*  *  olu— vlJ, 


for  every  u,v  €  V. 


(C)  (Equi-boundcdncss):  There  exist  a  constant  6  >  0,  independent  of 


q  €  Q  such  that 

IA(q)vl.  <  0(lvl  +  1), 
for  every  v  €  V. 

For  each  q  €  Q,  let  f(-;q)  €  LjfO.TiH)  and  u°(q)  €  H  and  assume  that 
the  mapping  q  -  u°(q)  is  continuous  from  Q  C  Q  into  H  and  that  the 
mapping  q  ~  f(t;q)  is  continuous  from  Q  C  Q  into  H  for  almost  every  t  € 
[0,T].  Also,  for  every  z  €  Z,  let  u  *•  $(u;z)  be  a  continuous  map  from 
C(0,T;H)  into  R+. 

We  consider  parameter  identification  or  inverse  problems  of  the  form: 

(ID)  Given  observations  z  €  Z,  determine  parameters  q  €  Q  which 
minimize 

ftq)  -  *(u0(q);z) 

where  u0(q)  «  u0(-;q)  is  a  mild  solution  to  the  initial  value  problem 

(3.1)  u(t)  +  A(q)u(t)  -  f(t;q),  0  <  t  <  T, 

(3.2)  u(0)  »  u°(q) 

corresponding  to  q  €  Q. 

By  a  mild  solution  to  (3.1),  (3.2)  we  mean  a  solution  in  the  sense  of 
Theorem  2.1.  To  be  more  precise,  for  each  q  €  Q  we  define  the  operator 
A0(q):  Dom(A0(q))  C  H  -»  H  to  be  the  restriction  of  the  operator  A(q)  to  the 
subset  of  V  given  by  Dom(A0(q))  -  (v  €  V:  A(q)v  €  H),  and  prove  the 
following  theorem. 


Theorem  3.1.  For  each  q  €  Q  the  operator  A0(q);  Dom(A0(q»  C  H  -•  H  is 
densely  defined  and  the  operator  Afl(q)  +  ul  is  m-accretive. 

Proof.  We  first  show  that  for  each  q  €  Q  the  operator  A(q)  +  ul:  V  -*  V*  is 
coercive.  If  {vj  C  V  with  lim  Iv  I  ■  •,  then  from  assumptions  (B)  and  (C)  we 

n  n-«  " 

obtain 

lim  <(A(q)  +  wl)vn,vn>/lvnl 
n*4® 

=  lira  {(<A(q)vn  -  A(q)0,vn>  +  u|vn|J)/lvnl  +  <A(q)0,vn>/lvnl) 

>  Hm  ("IV'I’/lv.l  -|<A(q)0,vB>|/lvnl) 

>  lim  {«lv  I  —  IA(q)0l.}  >  lim  «lv  I  —  B  -  “ 

n-**”  n  n-**  “ 

where  0  denotes  the  zero  vector  in  V.  It  follows  that  for  each  X  >  0,  the 
operator  I  +  x(A(q)  +  ul):  V  -*  V*  is  monotone,  everywhere  defined  on  V, 
hemicontinuous,  and  coercive.  Consequently  R(I  +  X(A(q)  +  ul))  -  V*  (see 
Barbu  [6],  Theorem  11.1.3)  and  therefore  X(I  +  X(A0(q)  +  ul))  -  H.  Also,  for 
u,v  €  Dom(A0(q)),  we  may  use  assumption  (B)  to  conclude 

£l  +  ~y  j|u-v|*  <  |u- v|*  +  Xalu-vl* 

<  |u-v|*  +  X<(A(q)  +  ul)u  -  (A(q)  +  ul)v,  u-v> 

-  <(I  +  V(A(q)  +  ul))u  -  (I  +  X(A(q)  +  ul))v,u-v> 

<  KI  ♦  UA0(q)  +  ul)u  -  (I  +  X(A„(q)  +  ul)v|  |u-v| 
or 

|u-v|  «  |u  -  V  +  X((A0(q)  +  ul)u  -  (A0(q)  +  ul)v)| 

which  proves  that  A0(q)  +  ul  is  m-accretive  on  Dom(A0(q))  c  H. 

To  show  Dom(  A0(q) )  -  H,  we  let  u  €H  and  for  each  n  -  1,2, ...  we  set 
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un  *  J(l/n;A0(q)+wI)u  €  Dom(A0(q)).  Then,  arguing  as  we  have  above,  we 
find 

|un|*  +  (I/n)aluBl*  <  <u  -  (l/n)A(q)0,un> 

<  |u|  |un|  +  (l/n)IA(q)fll,lunl 

where  8  is  again  the  zero  vector  in  V.  But  then 

(3.3)  d/2)|u/  +  (l/n)(a/2)lunIJ  <  (l/2)|u|l  +  (l/nKl/2a)IA(q)8l* 

<  (l/2)|u|J  +  (I/nKes/2a), 

from  which  it  immediately  follows  that  the  un  are  uniformly  bounded  in  H. 
Indeed,  from  (3.3)  we  see  that  (l/n)lunl*  and,  hence  luBl/,m,  is  bounded  so 
that  luBl/n-0  as  n-.** 

Also,  assumption  (C)  yields 

luB  -  ul.  -  (l/n)l(A0(q)  +  ul)unl.  <  (l/n){(B  +  U|t*)lunl  +  B). 

Since  the  last  term  in  the  estimate  above  tends  to  zero  as  n  -*  •,  we  find 
un  u  in  V*  as  n  -•  This,  together  with  the  fact  that  V  is  dense  in  H 
imply  that  uB  -*  u  weakly  in  H  as  n  -*  •  from  which  Dom(A0(q))»H 
immediately  follows. 

In  light  of  Theorem  3.1,  we  may  apply  Theorem  2.1  with  X  ■  H,  A  * 
A0(q)  and  f  -  f(  ;q).  We  conclude  that  there  exist  a  unique  nonlinear 
evolution  system  (U0(t,s;q):  0  <  s  <  t  <  T)  on  H  satisfying  (i),  (ii)  and  (iii). 
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The  mild  solution  u0(  -;q):  [0,T]  -*H  to  the  initial  value  problem  (3.1),  (3.2)  is  given 
by  u„(t;q)  -  Uo(t,0;q)u°(q)  for  t  €[0,T]. 

Remark.  Under  additional  hypotheses  on  f(-;q)  and  u°(q)  other  existence 
results  can  be  applied  to  obtain  somewhat  different  notions  of  a  solution  to 
the  initial  value  problem  (3.1),  (3.2).  For  example  (see  [6,  p.140-144])  if 
f(-;q)  €  W1,1(0,T;H)  and  u°(q)  €  Dom(A0(q)),  then  there  exists  a  unique  u(-,q): 
[0,T]  -  V  satisfying  u(  ;q)  €  W^O.T;!!),  A(q)u(  ;q)  €  L.(0,T;H)  and  u(t;q)  + 
A(q)u(t;q)  -  f(t;q)  a.e.  t  €  [0,T].  Or,  if  u°(q)  €  H  and  f(  ;q)  €  L,(0,T;V*) 
then  there  exists  a  unique  u(-;q)  which  is  V**valued  absolutely  continuous 
almost  everywhere  on  [0,T],  u(  ;q)  €  C(0,T;H)  n  L2(0,T;V),  u(  ;q)  €  L2(0,T;V*) 
and  u(t;q)  +  A(q)u(t;q)  -  f(t;q),  a.e.  t  €  (0,T).  If,  in  addition,  the  mapping  t  - 
t?f'(t;q)  is  an  element  in  L2(0,T;V*)  for  some  7*1,  then  the  mapping  t  -» 
t?u(t;q)  is  in  L2(0,T;V)  n  L«(0,T;H).  In  particular,  when  f(  *;q)  *  0,  the 
nonlinear  semigroup  (S„(t;q):  0  <  t  <  T)  on  H  defined  by  S0(t;q)  -  UQ(t;0;q), 
t  €  [0,T],  with  generator  -A0(q)  behaves  like  a  holomorphic  linear  semigroup 
in  that  it  smooths.  That  is,  S0(t;q)u°(q)  €  Dom(A0(q)),  t  €  (0,T],  and  the 
mapping  t  -*  t  S(t;q)u°(q)  is  an  element  in  LJ0,T;H)  for  every  u°(q)  €  H. 
Also,  some  generalizations  are  possible.  For  example,  in  assumption  (B),  the 
term  alu  -  vl*  can  be  replaced  by  a  term  of  the  form  a(lu-vl)lu  -  vl  where 
a(  •)  is  a  continuous,  strictly  increasing  function  on 
(0,“)  satisfying  «(0)  «  0  and  lim  a(x)  *  •.  Or,  the  terms  lu  -  vl1  in  (B) 

X-*” 

and  Ivl  in  (C)  can  be  replaced  by  lu  -  vlp  and  IvlP*1,  respectively,  for 
any  p  >  2. 
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The  development  of  computational  methods  for  the  solution  of  the 
infinite  dimensional  optimization  problem  (ID)  requires  the  finite  dimensional 
approximation  of  the  abstract  initial  value  problem  (3.1),  (3.2).  The  general 
framework  that  we  are  proposing  is  based  upon  a  classical  Galerkin  approach. 
For  each  n  -  1,2,  ...  let  Hn  denote  a  finite  dimensional  subspace  of  H  which  is 
a  subset  of  V.  Let  Pn:  H  -•  Hn  denote  the  orthogonal  projection  of  H  onto  Hn 
with  respect  to  the  <•,•>  inner  product.  We  assume  that  the  approximating 
subspaces  Hn,  and  the  projections  Pn  satisfy 

(D)  For  each  v  €  V,  limlP  v  -  vi  «  0. 

n-*»  n 

Note  that  assumption  (D)  and  V  densely  and  continuously  embedded  in  H 
imply  that  lim  |P  u  -  u|  ■  0  for  each  u  €  H. 

n-** 

For  each  q  €  Q  and  n  -  1,2,  ...  we  define  the  single-valued  operator  An(q): 
Hn  -  Hn  by  An(q)un  -  v„  for  un  €  Hn  where  vn  satisfies 

<A(q)un,wn>  -  <vn,wn>,  wn  €  Hn. 

That  An(q)  is  a  well  defined  operator  from  Hn  into  Hn  follows  from  the  Riesz 
Representation  Theorem  applied  to  the  Hilbert  space  Hn  and  the  bounded 
linear  functional  <A(q)un,->  on  Hn.  Also,  define  fn(-;q):  [0,TJ  -*  Hn  and 
u°(q)  €  Hn  by  fn(t;q)  -  P„f(t;q),  0  <  t  <  T,  and  u°(q)  -  Pnu°(q),  respectively. 
Note  that  f„(;q)  €  Ll(0,T;H|))  C  L^O.TjH)  and  that  |fn(t;q)|  <|f(t;q)|  for  q  €Q 
and  almost  every  t€  (0,TJ. 

We  consider  the  sequence  of  approximating  identification  problems  given 
by: 


(IDn)  Given  observations  z  €  Z,  determine  parameters  qn  €  Q  which 
minimize 
*n(q)  -  *(un(q);z) 

where  uD(q)  ■  un(  ;q)  is  a  mild  solution  to  the  initial  value  problem 
in  H 

n 

(3.4)  u„(t)  +  AB(q)un(t)  -  fn(t;q),  0  <  t  <  T 

(3.5)  un(0)  -  u°(q) 
corresponding  to  q  €  Q. 

From  the  definition  of  the  An(q)  and  the  assumptions  (B)  and  (C)  on  A(q), 
using  arguments  analogous  to  those  used  to  prove  Theorem  3.1,  it  can  be 
shown  that  the  operators  An(q)  +  ul  are  m-accretive  on  Hn.  It  then  follows 
from  Theorem  2.1  that  for  each  n  -  1,2,  ...  there  exists  a  unique  nonlinear 
evolution  system  (Un(t,s;q):  0  <  s  <  t  <  T}  on  Hn  satisfying  (i)  -  (iii)  in  the 
statement  of  that  theorem  with  X  -  Hn,  f(t)  «  fn(t;q),  and  =  uJ(q).  The 
mild  solution  to  the  initial  value  problem  (3.4),  (3.5)  is  given  by  un(t;q)  - 
Un(t,0;q)u“(q),  t  €  [0,T]. 

If  we  assume  for  the  moment  that  the  approximating  identification 
problems  (IDn)  have  solutions  qn  €  Q,  then  it  is  desirable  that  they  in  some 
sense  approximate  a  solution  q  to  the  original  identification  problem  (ID). 
This  is  in  fact  the  case.  For  suppose  that  it  can  be  shown  that  for  any 
sequence  (qn)  C  Q  with  lim  qn  =  q  €  Q  we  have 

(3.6)  j»mun(qn)  -  u0(q0)  in  C(0,T;H). 

Then  (qn)  C  Q  and  Q  a  compact  subset  of  the  metric  space  Q  imply  that  there 


exist  a  subsequence  (qn }  C  (qn)  and  a  q  €  Q  such  that  lim  qn  -  q.  For  any 
q  €  Q  the  continuity  of  *  implies 

#q)  -  4<u0(q);z)  -  *  ( lim  u  (q  );z) 

0  j—  “j  “j 

-  lim  *(u  (q  );z)  -  lim  4>  (q  ) 

J-**  "j  n j  j—  "j  Oj 

<  lim  6  (q)  -  lim  6(u  (q);z) 
i"**  j  j-*  j 

-  *  (lim  u  (q);z)  -  *(u0(q);z) 

j"**  j 

-  «(q). 

Note  that  in  the  discussion  above  we  did  not  assume  that  a  solution  to 
problem  (ID)  exists.  But  rather  we  have  shown  that  the  existence  of  solutions 
qn  to  the  approximating  problems  (IDa)  and  (3.6)  imply  the  existence  of  a 
solution  q  to  problem  (ID).  When  the  solution  to  problem  (ID)  is  unique,  the 
sequence  (qn)  itself  converges  to  q. 

The  existence  of  a  solution  qB  to  problem  (IDn)  for  each  n  -  1,2,  ...  will 
follow  from  the  compactness  of  Q  and  the  continuity  of  *  once  the  continuous 
dependence  result:  lim  un(qm)  -  u„(Qo^  ‘n  C(0,T;Hn)  whenever  (qm)  C  Q  with 
lim,  qm  *=  qQ,  has  been  established.  Although  continuous  dependence  for  the 
finite  dimensional  systems  (3.4),  (3.5)  could  be  demonstrated  via  a 
modification  to  any  one  of  a  number  of  familiar  continuous  dependence 
results  for  ordinary  differential  equations  (see,  for  example.  Hale  [12], 
Theorem  1.3.4),  it  is  also  easily  handled  with  the  approximation  theory 
developed  in  the  previous  section.  This  and  the  convergence  in  (3.6)  are 
addressed  in  the  following  theorem. 
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Theorem  3.2.  If  assumptions  (A)  -  (D)  hold ,  then 

(a)  If  {q  }  C  Q  with  lim  qn  -  q0  then  lim  u  (q  )  =  u0(q0)  in  C(0,T;H),  and 

n“*®  n“*® 

(b)  If  (qm}  C  Q  with  lim  qm  -  q0  then  lim  un(qm)  =  un(q0)  for  each  n  €  Z+. 

m~*°°  m“*°° 

Proof.  Assumption  (D)  and  the  continuity  of  the  map  q  -»  u°(q)  from  Q  into  H 
imply  lim  u°(qn)  -  u°(q0)  in  H.  Hence,  we  will  have  verified  (a)  if  we  can 
show  that  lim  Un(t,s;qn)wn  »  U0(t,s;q0)wQ,  0  <  s  <  t  <  T,  uniformly  in  t  for 

n"*® 

t  €  [s,T]  whenever  wn  €  Hn  with  lim  wn  »  wQ  €  H.  We  argue  this  using 

n-*® 

Theorem  2.2.  Note  that  assumption  (D)  implies  lim  H  D  H  and  assumption 

n-**  " 

(D)  together  with  the  assumed  continuity  of  the  map  q  -*  f(t;q)  from  QC0 
into  H  for  almost  every  t  €[0  ,T]  imply  lim  f  (t  ;q  )  ■  f(t;qn)  in  H  for  almost 
every  t€  [0,T]  with  the  fn(';qn)  dominated  by  a  function  g€Lj(0,T;H) 
which  is  independent  of  n.  Thus,  we  need  only  to  demonstrate 
that  for  some  XQ  >  0  we  have 

(3.7)  limJ(X0;An(qn)  +  wl)wn  =  JfX^A^q,,)  +  ul)w0 


in  H  whenever  w  €  H  ,  n  €  Z+  with  lim  w  =  wn. 

n  n  n  u 


Let  Xfl  >  0  and  set  vn  -  J(X0;An(qn)  +  ul)wn  and  vQ  -  J(X0;A0(q)  +  wl)w0. 
We  first  show  that  Iv  I  is  uniformly  bounded  in  n.  From  assumption  (B)  we 


obtain 


X0alvn|1  <  +  X0<A^n)V»  “  A(‘l„)0’V„> 

-  <(I  +  X0(An(qn)  +  WI))vn,vn>-|v/ 

+  x0 <A(QO)0  -  A«Oe.V„>  -  XO<A^o)0*Vn> 

“  <w»«vn>  "  lvJ*  +  X0<A(q0)®  "  A(£«n)0’Vn>  “  XO<A<£»O)0’Vn  > 

<  lwnl.lvnl  +  XoiA(qo)0  -  A(qn)0Mvnl  +  XolA(qo)0l.lvnl 
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where  0  denotes  the  zero  vector  in  V.  This  estimate  together  with  assumption 
(C)  yields 

lvnl  <  (X0«)'1K|wn|  +  «_1||A(qn)0  -  A(qo)0l.  +  a_16. 

Recalling  assumption  (A)  and  that  lim  w  =  wrt  in  H,  we  find  that  the  desired 

n-**  n  0 

uniform  bound  on  Iv  I  has  been  established. 

n 

Once  again,  from  assumption  (B),  we  find 

X0alv„^0|J  «  +  X0<A(q„)V„  "  A(^n)V0*Vn  "  V0> 

=  V0UlVn-Vol2  +  V^nK  "  A(%>V0’Vn  "  P„V0> 

+  X0<A(<UVn  "  A(£>0)V0*PnV0  "  V 
+  Xo<A^o)V0  “  A<<ln)v0,Vn  ~  V0> 

"  X0U<PnV0  ~  V0’V„  "  V0>  +  «!  +  Xo(An<(ln)  +  Ul»V„ 

-  (I  +  yA0(q0)  +  wl))v0,vn  -  PnvQ>  +  <vQ  -vn,vn-Pnv0> 

+  Xo<A^n)Vn  ~  A<<»0>vO-PnV0  "  V0> 

+  X0<A<<»o)V0  “  A<<»n>vO*Vn  ~  V0> 

=  X0U<PnV0  _  V0’Vn  “  V  +  <Wn  _  W0’Vn  ~  PnV0>  “I  v„  _P„V0  I* 

+  X0<A(<ln>Vr.  ~  A^o)VPnV0  “  V 

+  X0<A^O>VO  ~  A^n>v0*vn  "  V0> 

<  X0^IPnV0_ V0I* 1  Vn~ V0^  +  ,w„-w0,‘lv^0l 

+  lwn-w0IJPnv0^0l  +  X0lA(qn)vn  -  A(q0)v0I.IPnv0  -  vQl 

+  X0|A(%)V0  "  A(flJV0,<,Vn  _  Vo'- 
1  .  n  , 

The  estimate  ab  <  —  a*  +  -  b  for  any  n  >  0  and  assumption  (C)  allow  us  to 
2n  2 


argue 


~2~’ Vn~”V0|2  *  ^T^Vo^O1*  +  ^77  ,Wn~ W0I*  +  l*n^0l'IPnV0^0l 


+  X0,A<<ln)Vn  ~  A(<lo)V0,*,PnV0  “  Vo'  +  ^f-1  A(%)V0  "  A^K1* 


3w2\0m4  (X0+l) 


}  'Vo  "  Vo‘’  +  {  *  Y  )|Wn  "  Wol1 


+  \,{«lvnl  +  lv0*)  +  2>,PnV0  "  Vo'  +  ,A(Qo)V0  “  A(qn)Vo'*  • 


From  this,  the  uniform  bound  on  lvnl,  lim  wn  =  wQ  in  H  and  assumptions 
(A)  and  (D)  allow  us  to  conclude  lim  vn  =  vQ  in  V  and  that  (3.7)  holds. 

An  analogous,  but  somewhat  simpler  argument  can  be  used  to  verify  (b). 
We  use  Theorem  2.2  to  show  that  for  n  €  Z+  fixed,  lim  U  (t,s;q  )w  = 

—  n  m  m 


Un(t,s;q0)w0,  0  <  s  <  t  <  T,  uniformly  in  t  for  t  €  [s,T]  whenever  wm,wQ  €  Hn 
with  -  w0  in  H  Clearly  Un^f^^qJ  =  fn(  ;q0)  in  Lx(0,T;Hn)  so  that 

we  need  only  to  show  that  for  some  X0  >  0, 

+  Ul>Wm  "  +  Ul>W0 

in  H  whenever  Hm  wm  -  wQ  in  H.  Let  vm  -  J(X0;An(qJ  +  wl)wm  and  vQ  = 
J(  xo»An^qo^  +  W*)w0.  Then  from  assumption  (B) 

x0a,v«“v0l*  <  X0UlVm“Voll  +  X0<A(qm)Vm  ~  A(qm>  VVm”  V0> 

-  <(I  ♦  X0(An(qm)  +  ul))vro  -  (I  +  X0(An(q0)  +  U)I))v0,vm  -  v0> 

"  lvnTvJ*  +  xo<A(qo)vo  “  A(qm>Vvm  “  vo> 

"  <WnT"v  Vm“V0>  “  +  X0<A(q0>V0  “  A(qm>v0’Vm  “  V 

<  1  1  * 1  Vm~V0 1  +  X0,A<q0>V0  "  AK„>V0,*,Vm  "  V 


IVm  "  V0I  <  77  lWm  -  W0>  +  l~  'A(VV0  "  AfoJV* 
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Assumption  (A)  and  lim  w  =  w.  in  H  yield  the  desired  result  and  the 

m-*  m  0 

theorem  is  proved. 

Remark.  In  practice,  the  approximating  identification  problems  (IDn)  are  solved 
using  standard  iterative  search  techniques  (for  example,  steepest  descent,  Newton’s 
method,  etc.)  requiring  the  evaluation  of  4>n(q)  for  q  €Q  at  each  step.  This  in  turn 
requires  the  integration  of  the  finite  dimensional  initial  value  problem  (3.4),  (3.5). 
Once  a  basis  for  Hn  has  been  chosen,  the  solution  to  (3.4),(3.5)can  be  computed  using 
any  standard  numerical  integrator  for  ordinary  differential  systems.  Also,  the 
parameter  space  Q  and  the  admissible  parameter  set  Q  are  frequently  functional  in 
nature  and  infinite  dimensional.  When  this  is  the  case  the  set  Q  must  also  be 
discretized.  Suppose  that  for  each  m  =  1,2,...,  1“:  Q  CQ  -Q  is  a  continuous  map  with 
finite  dimensional  range  and  that  lim  I^q)  *  q  with  the  convergence  uniform 

m"*®“ 

in  q  for  q  €  Q.  Set  Qm  =  I^Q)  (note  that  Qm  is  a  compact  subset  of  Q)  and 
consider  the  identification  problems  (I D®)  defined  to  be  the  problems  (IDn) 

with  Q  replaced  by  Qm.  It  is  clear  that  each  of  these  problems  admit  a 

solution  q^1  and  it  is  not  difficult  to  argue  that  there  exists  a  subsequence 

{q^}  C  (q^1)  with  lim  q^j-  q,  q  a  solution  to  problem  (ID)  (see,  for  example, 
k  j.k-**  "k 

[4]).  Once  bases  for  Hn  and  the  range  of  Im  have  been  chosen,  problem  (ID™) 
involves  the  minimization  of  a  functional  over  a  compact  subset  of  Euclidean 
space  subject  to  finite  dimensional  constraints. 

Remark  (Nonautonomous  systems).  Theorems  2.1  and  2.2  remain  valid  for 
certain  classes  of  temporally  inhomogeneous  or  time  dependent  operators  A  = 
A(t).  To  be  more  precise,  the  family  of  operators  A(t):  X  -*  2X  must  be 
m-accretive  on  X  for  almost  every  t  €  [0,T]  and  must  satisfy 


(3.8)  |J(X;A(t))x  -  J(X;A(s))x|x  <  X|h(t)  -  h(s)|xL(|x|x) 


for  each  x  €  X,  every  X  satisfying  0  <  X  <  XQ  for  some  XQ  >  0,  some  h  € 
Lj(0,T;X),  some  continuous,  non  decreasing  function  L:  [0,«)  -  [0,*)  and 
almost  every  t,s  €  [0,T]  (see  [8],  [9]).  (Note  that  for  simplicity  we  have  taken 
u  =  0;  however,  the  discussion  to  follow  remains  valid  for  any  w  €  R.)  The 
primary  motivation  for  developing  the  framework  outlined  above  was  to 
define  readily  verifiable  conditions  on  the  operators  A(q):  V  -♦  V*  that  if 
satisfied  would  (i)  also  automatically  be  satisfied  by  the  Galerkin 
approximation  An(q)  and  (ii)  lead  to  the  desired  convergence  of  solutions  to 
the  approximating  identification  problems  to  a  solution  to  problem  (ID).  The 
natural  assumption  to  add  to  (A)  -  (C)  that  certainly  satisfies  criterion  (i)  and 
that  could  conceivably  lead  to  an  estimate  of  the  form  (3.8)  in  H  is  that 

(3.9)  lA(t;q)v  -  A(s;q)vl,  <  |h(t)  -  h(s)|L(|v|) 

for  each  v  €  V,  almost  every  s,t  €  [0,T]  and  some  h  €  Lj(0,T;H)  and  some 
continuous  nondecreasing  L:  (0,-)  -  {0,-),  both  of  which  do  not  depend  upon 
q  €  Q.  Unfortunately,  however,  we  can  only  show  that  (3.9)  leads  to  an 
estimate  of  the  form 

(3.10)  |J(X;A0(t;q))u  -  J(X;A0(s;q))u|  <  A  |h(t)  -  h(s)|  L(|u|) 

for  each  u  €  H.  Moreover,  it  is  not  clear  to  us  how,  or  if,  the  proof  of  the 
fundamental  Theorem  2.1  given  in  {9J  could  be  modified  so  that  (3.10)  would 
suffice.  We  have  explored  alternative  approaches  and  developed  other 
techniques  for  treating  the  nonautonomous  case  (for  example,  in  the  linear 
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case,  based  upon  some  ideas  in  Tanabe  [18],  and  in  the  strongly  monotone  case, 
via  a  variational  formulation  which  can  be  found  in  Barbu  (6]>.  These 
results  will  appear  soon  in  forthcoming  papers. 


4.  Applications  and  Examples 

We  briefly  describe  some  classes  of  systems  to  which  the  general 
framework  developed  in  the  previous  section  applies.  In  our  discussion  below 
we  consider  theoretical  aspects  only.  Implementation  questions  will  be  treated 
and  the  results  of  our  numerical  studies  will  be  reported  on  elsewhere. 

Example  4.1.  Linear  regularly  dissipative  operators.  The  approximation 
theory  for  inverse  problems  for  systems  involving  linear  regularly  dissipative 
operators  was  treated  in  detail  by  Banks  and  Ito  in,  and  is  the  central  focus 
of,  [2]  and  [3].  We  show  here  that  the  linear  theory  is  a  special  case  of  the 
nonlinear  theory  given  in  Section  3. 

Let  the  spaces  H,  V,  V*  and  Q  and  the  set  Q  be  as  they  have  been  defined 
above.  For  each  q  €  Q  let  a(q)(-,-)  be  a  sesquilinear  form  defined  on  V  x  V 
which  satisfies  the  conditions: 

(A1)  For  each  v  €  V  the  mapping  q  -*  a(q)(-,v)  is  continuous  from  Q  c  Q 
into  V*.  That  is  given  «  >  0  there  exists  a  6  >  0  such  that 

suj5  |a(qQ)(u,v)  -  a(q)(u,v)|  <  t 

lul=l 

whenever  d(q0,q)  <  6  where  d  denotes  the  metric  on  Q. 

(B')  There  exist  an  w  €  R  and  an  o  >  0,  both  independent  of  q  €  Q,  for 
which  a(q)(v,v)  +  ujv|*  >  olvl2  for  every  v  €  V. 

(C')  There  exists  a  constant  8  >  0,  independent  of  q  €  Q,  such  that 

|a(q)(u,v)|  <  Blul  Ivl  for  every  u,v  €  V. 


When  conditions  (A1)  -  (C*)  are  satisfied  it  is  not  difficult  to  argue  that  for 
each  q  €  Q  an  operator  A(q)  €  t(V,V*)  can  be  defined  by 

[A(q)v](u)  -  <A(q)v,u  >  -  a(q)(u,v), 

u,v  €  V  and  that  A(q):  V  -»  V*  satisfies  (A)  -  (C).  It  then  follows  from 
Theorems  2.1  and  3.1  that  there  exists  a  unique  nonlinear  evolution  system 
(U0(t,s;q):  0  (  s  <  t  <  T)  on  H  corresponding  to  the  initial  value  problem 

u(t)  +  A0(q)u(t)  -  f(t;q),  0  <  t  <  T 
u(0)  -  u°(q) 

where  for  each  q  €  Q,  f(;q)  €  1^(0, T;H),  u°(q)  €  H  and  A0(q):  Dom(A0(q))  C 
H  -*  H  is  the  restriction  of  A(q)  to  the  set  Dom(A0(q»  «  {v  €  V:  A(q)v  €  H}. 
The  operator  -AQ(q)  is  the  infinitesimal  generator  of  an  analytic  semigroup 
(T0(t;q):  t  >  0}  on  H  (see  [18])  and  for  $  €  H 

(4.1)  U0(t,s;q)*  -  T0(t-s;q)*  +  [  T0(t-r;q)f(T;q)dT. 

•i 

It  can  be  shown  that  the  semigroup  (T0(t;q):  t  >  0}  admits  an  extension  (T(t;q): 
t  >  0)  which  is  an  analytic  semigroup  on  V*  with  generator  A(q):  V  c  V*  -*  V*. 
Also  the  restriction  of  (T0(t;q):  t  >  0)  to  V,  call  it  (T(t;q):  t  >  0),  is  an  analytic 

A#  A 

semigroup  on  V  with  generator  A(q):  Dom(A(q))  CV-*V,  the  restriction  of 
A(q)  to  the  set  Dom(A(q))  ■  {v  €  V:  A(q)v  e  V)  (see  {3],  [18]).  Consequently, 
with  appropriate  assumptions  on  f(-;q),  the  evolution  system  (U0(t,s;q):  0  <  s 
<  t  <  T)  admits  an  extension  (U(t,s;q):  0  <  s  <  t  <  T)  which  is  an  evolution 

A* 

system  on  V*  and  a  restriction  (U(t,s;q):  0  <  s  <  t  <  T)  which  is  an  evolution 
system  on  V. 
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It  is  clear  from  (4.1)  that  when  A(q)  is  linear,  we  may  take  f(-;q)  =  0  and 
consider  only  the  approximation  of  the  semigroup  (T0(t;q):  t  *  0).  For  each 
n  -  1,2,  ...  let  the  finite  dimensional  subspaces  Hn  of  H  and  the  corresponding 
orthogonal  projections  Pn  be  as  they  were  defined  in  Section  3  and  assume 
that  condition  (D)  is  satisfied.  Denote  the  Galerkin  approximations  to  A(q) 
(i.e.  the  restriction  of  A(q)  to  an  operator  from  Hn  into  H*  =  Hn)  by  An(q)  and 
set  Tn(t;q)  -  exp(-tAn(q)),  t  *  0.  Theorem  3.2  then  implies  that 

(4.2)  lim  |Tn(t;qn)Pnu°(qn)  -  T0(t;q0)u°(q0)|  -  0 

uniformly  in  t,  for  t  €  [0,T]  whenever  {q  }  C  Q  with  lim  q  -  q0  €  Q,  and 
the  mapping  q  -  u°(q)  is  continuous  from  Q  C  Q  into  H.  In  addition,  recalling 

that  we  required  that  Hn  C  V  for  all  n  -  1,2 . an  inspection  of  the  proof  of 

Theorem  3.2  reveals  that  in  the  linear  case  with  the  existence  of  the  semigroup 
(T(t;q):  t  >  0}  on  V,  we  may  apply  Theorem  2.2  with  X  «  V  and  conclude  that 

(4.3)  Urn  lTn(t;qn)Pnu°(qn)  -  T  (t;q0)u°(q0) I  -  0 

uniformly  in  t  for  t  €  [0,T]  whenever  lim  q  -  q^  u°(q)  €  V  and  the  map 

it4* 

q  -•  u°(q)  is  continuous  from  Q  into  V  (see  also  [3]).  Then  for  $  €  H,  setting 

Un(t,s;q)Pn*  *  Tn(t-s;q)Pn*  +  f  Tn(t-r;q)PRf(T;q)dT 

*  • 

under  appropriate  assumptions  on  f(-;q),  (4.2)  and  (4.3)  continue  to  hold  with 
I  Tn(t;q),  T0(t;q),  and  T(t;q)  replaced  by  UB(t,s;q),  U0(t,s;q),  and  U(t,s;q), 

I  respectively,  with  the  convergence  being  uniform  in  t,  for  t  €  (s,TJ.  Hence  the 

i 

I 

i  linear  theory  and  results  of  (3]  are  a  special  case  of  the  nonlinear  theory  of 

t 


Section  3. 
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We  note  that  in  the  context  of  the  identification  problem,  the  fact  that 

the  stronger  V-convergence  given  in  (4.3)  can  be  obtained  is  significant. 

Indeed,  (4.3)  permits  the  relaxation  of  the  continuity  assumption  on  the 

performance  index  ♦  to  the  requirement  that  for  each  z  €  Z,  the  mapping  u  - 

4<u,z)  be  continuous  from  C((0,T];V)  into  R+.  This  can  have  the  effect  of 

significantly  enlarging  the  class  of  allowable  observations.  For  example,  in 

the  case  of  a  one  dimensional  parabolic  system  formulated  in  H  «  L2  with  V 

in  Hl,  spatially  discrete  (i.e.  pointwise,  as  opposed  to  distributed  in  space) 

measurements  will  suffice  (see  [3]  and  [3]). 

Among  the  class  of  linear  regularly  dissipative  operators  which  arise  from 

a  form  satisfying  (A1)  •  (C1)  are  the  familiar  elliptic  partial  differential 

•  ls+f+i 

operators  on  Lr  Briefly,  let  A  be  a  region  in  R'  and  let  Q  -  *  L„(A).  Let 

n=l 

Q  be  a  compact  subset  of  Q  with  the  property  that  if  q-  {(ay,  bs,  c):  i,j«J,...,f]  €  Q, 
then  for  some  a  >  0  independent  of  q  €  Q. 


for  every  x  €  A,  and  every  (  €  R*.  For  q  €  Q  and  u,v  €  H*(A)  set 


a(q)(u 


•v>-  UL ‘“<x) 


du(x)  dv(x) 


dx.  dx 


i 


1  3u(x)  > 

+  E  b.(x)  - - -  v(x)  +  c(x)u(x)v(x)}dx 

i=i  1  ox.  J 


I 


with  H  -  L2(A)  and  V  any  closed  subspace  of  HHA)  containing  h£(A),  it  can  be 
shown  (see  [18])  that  a(q)(-,  )  satisfies  (A1)  -  (C').  The  operator  A(q)  is  given 
formally  by 
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►  *«*  1.^  4L  «.'  *» 


(4.4)  A(q)  -  -  E  aij(x)  I-  +  .£  Vx>  I- +  c(x>- 

ij=i  9Xj  flXj  i=i  1  8Xj 

When  do  is  sufficiently  smooth,  A(q)  is  the  elliptic  operator  given  by  (4.4),  and 
V  is  chosen  to  be  either  Hj(ft)  or  Hx(ft),  the  equation  (3.1)  becomes  a  parabolic 
partial  differential  equation  with  either  Dirichlet  or  Neumann  boundary 
conditions. 

For  H  ■  Lj(fl)  and  V  a  subspace  of  H*(ft),  choosing  the  approximating 
subspaces  to  be  the  span  of  an  appropriate  collection  of  first  order  spline 
functions  will  typically  satisfy  assumption  (D)  (see  [IS]  and  Example  4.2 
below). 

Example  4.2.  Nonlinear  Elliptic  Operators.  Let  ft  be  a  bounded  region  in  R* 
with  smooth  boundary  r  «  dft.  For  a  *  («r  at)  a  multi-index,  let  |a|  *  + 
<*j  +  •  •  •  +  <*|  and  denote  the  ath  order  generalized,  or  distributional 
derivative  of  a  function  u  by  D®u;  that  is. 


D“u(x) 


— —  u(x),  x  €  ft  . 

V 


Let  m  be  a  nonnegative  integer  and  let  Su  denote  the  vector  valued  function 
of  length  N  -  (,Jm)  whose  components  are  all  of  the  partial  derivatives  of  u 
of  order  greater  than  or  equal  to  zero  and  less  than  or  equal  to  m. 

For  each  multi-index  a  with  M  <  m,  let  (x,S)  -»  a^x.C)  be  a  real  valued 
function  defined  on  ft  *  RN  which  is  measurable  in  x  and  continuous  in  We 
assume  that 

(1)  there  exist  a  g  €  L2(ft)  and  a  positive  constant  y  such  that 


V 1 »7r\# V ‘Tc1  » m 


yv.m  M 
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(4.5)  la^x.OI  <  7(ltl  +  g(x)) 

for  almost  every  x  €  Cl,  each  (  e  RN  and  all  a  with  |<4  <  m,  and 
(2)  there  exists  a  positive  constant  X  such  that 


H'6)  »$.  k  wr«»  ll«  '  nJ’ 

for  almost  every  x  €  ft  and  all  {,n  €  RN. 

Let  H  =  L2(ft)  and  let  V  be  any  closed  subspace  of  H^ft)  which 
contains  H™(ft).  Define  the  operator  A:  V  -»  V*  by 


(4.7)  (Au)(v)  =  E  a^x,6u(x))D“v(x)dx, 


for  u,v  €  V.  The  operator  A  given  by  (4.7)  is  the  distributional  form  of  the 
formal  differential  operator 


(4.8)  (Au)(x)  =  E  (-l)aD«aa(x,8u(x)). 

|oq*m 

A  differential  operator  of  the  form  (4.8)  is  referred  to  as  a  nonlinear  elliptic 
operator  and  the  partial  differential  equation 

du 

(4.9)  —  (t,x)+  E  (~I)0fD°aa(x,8u(t,x))  -  f(t,x) 

Ot  |a|  <m 

is  said  to  be  of  nonlinear  parabolic  type.  When  V  -  H£\ft),  a  solution  in  V* 
to  the  abstract  equation 


u(t)  +  Au(t)  -  f(t) 

with  A  given  by  (4.7)  corresponds  to  a  variational  solution  to  (4.9)  which 
satisfies  Dirichlet  boundary  conditions.  When  V  »  H^ft),  a  variational 
solution  to  the  Neumann  problem  is  obtained.  Note  that  in  the  linear  case  we 
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have 

aJx.SuU))  -  E  a°^(*)D®u(x). 

|p|*m 

Under  the  assumptions  above,  it  is  not  difficult  to  show  that  A  given  by 
(4.7)  is  hemicontinuous  and  satisfies  conditions  (B)  and  (C)  given  in  Section  3. 
With  an  appropriate  choice  of  the  space  Q  and  the  set  Q,  condition  (A)  can  be 
satisfied  as  well. 

A  quasilinear  model  for  heat  conduction  or  mass  transfer  in  which  the 
heat  or  mass  flux  is  a  function  of  the  temperature  or  mass  fraction  gradient 
discussed  in  [16]  and  [17]  leads  to  a  nonlinear  elliptic  operator  and  a  nonlinear 
parabolic  partial  differential  equation  of  the  forms  (4.8)  and  (4.9),  respectively 
with  m  -  1.  Let  0  be  a  bounded  region  in  R1  with  smooth  boundary  and  let 
Q  «  Ljofft  x  R1).  Let  Q  be  a  compact  subset  of  Q  with  the  property  that 
q  €  Q  if  and  only  if  the  mapping  (  -*  q(x,t)  is  Cl  on  R1  for  almost  every 
x  €  fl  and  there  exists  a  X  >  0  (which  does  not  depend  on  q)  such  that 

(4.10)  ei^q(x,0|t=e  (t  -  0)  +  q(x,8Xt,  -  hj)  >  XUt  -  I).), 

for  i  -  1,2,  ....  *,  almost  every  x  €  n  and  all  0,t,>7  €  R1.  (When  1-1,  the 
function  q(x,l)  -  q(?)  -  (1  -  .5e'**)  satisfies  (4.10).) 

Let  H  -  Lj(O)  and  let  V  be  any  closed  subspace  of  Hl(0)  which  contains 
Hj(G).  Then  V  c  H  C  V*  and  for  each  q  €  Q  define  A(q):  V  -*  V*  by 

(4.11)  (A(q)uXv)  -  [  q(x,vu(x))Vu(x)- W(x)dx 

for  u,v  €  V.  Note  that  for  each  q  €  Q  the  operator  given  by  (4.11)  is  of  the 
form  (4.7)  with 
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(4.12)  a<x(x, 6u(x))  -  q(x,Vu(x))D“u(x) 

for  x  €  fl  and  all  «  with  M  -  1  and  aa  *  0  for  |a)  -  0.  The  nonlinear 
parabolic  partial  differential  equation  (4.9)  takes  the  form 

du 

—  (t,x)  -  7  •  q(x,Vu(t,x))VU(t,x)  «  f(t,x),  t  >  0,  x  €  n  . 

Taking  I  •  I  to  be  the  usual  norm  on  H^ft),  it  follows  that 
lA(q0)u  -  A(q1)ul*  <  |q0  -  q^lul 

for  each  u  €  V  and  q^!  €  Q.  Since  Q  is  a  compact  subset  of  x  R4),  it  is 
easily  verified  that  aa  given  by  (4.12)  satisfies  a  growth  condition  of  the  form 
(4.5)  with  7  and  g  independent  of  q  €  Q.  An  application  of  the  mean  value 
theorem  together  with  assumption  (4.10)  imply  the  existence  of  a  X  >  0, 
independent  of  q  €  Q,  for  which  (4.6)  holds.  Consequently  the  conditions  (A),  (B), 
and  (C)  given  in  Section  3  are  satisfied,  and  our  general  theory  can  be  applied. 

With  regard  to  approximation,  polynomial  spline  function  based  Galerkin 
subspaces  can  often  be  shown  to  satisfy  condition  (D).  For  example,  when 
I  •«  1  and  O  *  (0,1)  in  the  nonlinear  heat  conduction/mass  transfer  example 
discussed  above,  the  subspaces  Hn  can  be  chosen  as  the  span  of  the  linear 
B-spline  ("hat")  functions  with  respect  to  the  uniform  mesh  {0,  1/n,  2/n,  ...,  1} 
appropriately  modified  to  satisfy  stable,  or  geometric,  boundary  conditions. 
Familiar  error  estimates  for  interpolation  and  the  Schmidt  inequality  can  then 
be  used  to  verify  that  condition  (D)  is  satisfied  (see  [5]).  Generalization  to 
higher  dimensions  is  possible,  and  can  often  be  achieved  via  tensor 
products  of  one  dimensional  elements  (see  [15]). 


5.  Concluding  Remarks 


We  have  developed  a  general  abstract  approximation  framework  for  the 
identification  of  nonlinear  distributed  parameter  evolution  systems.  The  class 
of  systems  to  which  our  theory  applies  are  those  whose  dynamics  can  be 
described  by  a  nonlinear  operator  which  satisfies  conditions  that  are  the 
natural  nonlinear  extensions,  or  analogs,  of  the  properties  of  regularly 
dissipative,  or  abstract  parabolic,  linear  operators.  The  approach  we  have 
taken  is  based  upon  the  defining  of  a  sequence  of  approximating  finite 
dimensional  identification  problems  in  which  the  systems  to  be  identified  are 
Galerkin  approximations  to  the  original,  underlying,  infinite  dimensional 
nonlinear  dynamics.  Under  a  weak  continuity  assumption  with  respect  to  the 
unknown  parameters  to  be  identified,  equi-boundedness  and  equi- 
monotonicity  conditions,  and  an  approximation  assumption  on  the  Galerkin 
subspaces  (all  of  which  are  readily  verified  for  wide  classes  of  nonlinear 
distributed  systems  and  finite  element  subspaces),  we  are  able  to  demonstrate 
that  solutions  to  the  approximating  problems  exist,  and,  in  some  sense, 
approximate  (i.e.  subsequential  convergence)  solutions  to  the  original  infinite 
dimensional  identification  problem.  We  have  shown  that  the  linear  theory 
presented  in  [2]  and  [3]  is  a  special  case  of  our  nonlinear  framework  and  that 
our  results  are  applicable  to  a  reasonably  wide  class  of  nonlinear  elliptic 
operators  and  corresponding  nonlinear  parabolic  partial  differential  equations. 
In  particular,  we  have  considered  application  of  our  theoretical  framework  to 
a  quasi-linear  model  for  heat  conduction  or  mass  transport. 


The  general  approximation  result  for  nonlinear  evolution  systems 
discussed  in  Section  2  is  applicable  to  a  much  broader  class  of  nonlinear 
dynamical  systems  than  we  subsequently  treated  in  Section  3.  For  example, 
this  class  of  systems  would  include  those  with  dynamics  described  by  set 
valued  maps  or  multifunctions,  and  (after  minor  modification  to  the 
general  theory)  time  dependent  or  nonautonomous  operators.  We  are 
currently  investigating  these  features  of  the  general  approximation  theory 
in  the  context  of  parameter  estimation  problems.  Also,  we  would  like  to 
be  able  to  weaken  the  somewhat  restrictive  strong  monotonicity  condition. 
Any  progress  that  we  might  make  in  these  efforts  would  have  the  potential 
to  significantly  enlarge  the  class  of  nonlinear  systems  to  which  our  theory 
and  framework  would  apply.  Finally,  extensive  numerical  or 
computational  studies  designed  to  demonstrate  the  feasibility  and  point  out 
the  limitations  of  our  schemes  and  general  approach  are  currently 
underway  and  will  be  reported  on  in  a  forthcoming  paper. 
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